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1. INTRODUCTION 

The HERA collider circulates beams of elec- 
trons (positrons) and protons at energies of 27.5 
and 820 GeV respectively. The e^p centre of 
mass energy for collisions between the beams is 
300 GeV. The largest cross-section processes ob- 
served in these collisions are photoproduction re- 
actions, where the beam acts as a somxe 
of almost real {Q^ <C IGeV^) photons. Pro- 
cesses where a virtual photon of sizeable is 
exchanged are known as Deep Inelastic Scatter- 
ing, or DIS reactions. In this paper I review mea- 
surements of open charm production at HERA 
made in both photoproduction and DIS by the 
HI and ZEUS collaborations. 



a) Photon-gluon fusion PGF b) Charm flavour excitation 




c) Resolved photon contributions 



Figure 1. Charm production processes at 
HERA. 



*on behalf of the ZEUS and HI collaborations. 



Charm quarks are produced in real or virtual 
photon-proton collisions predominantly via the 
photon-gluon fusion process illustrated in Fig- 
ure Charm production is therefore sensitive 
to the gluon momentum distribution in the pro- 
ton. Other contributing processes are also shown 
in Figure |l]. The dia gram of Figure |l|b, where the 
charm quark is an active constituent of the pro- 
ton, is expected to be important in DIS at large 
Q^. In photoproduction, additional contributions 
may arise from "resolved photon" processes as 
shown in Figure [^c, which reveal the partonic 
structure of the photon. Sources of charm pro- 
duction other than via the hard processes shown, 
such as fragmentation or the decay of b quarks, 
are negligible by comparison. 

Cross-section calculations for these processes in 
the framework of perturbative QCD should be re- 
liable, since the relatively large charm quark mass 
rric > Aqcd provides a hard scale in all cases. 
In addition, in many kinematic regions, further 
hard scales are provided by a large photon vir- 
tuality or high transverse momentum of the 
produced charm quark pair. A number of authors 
have produced calculations of charm production 
to next-to-leading order (NLO) in QCD. 

The analyses reviewed here are based on data 
taken in 1994, when each experiment collected 
around 3 of positron-proton collisions. Both 
HI and ZEUS reconstruct charm signals using the 
invariant masses of tracks found in the tracking 
detectors. No particle identification is used. I 
cover only the production of open charm; cc me- 
son signals are discussed in the contribution to 
this conference on vector mesons at HERA[0. 

Throughout this paper I use standard symbols 
for event-related kinematic quantities at HERA, 
is the invariant squared 4-momentum transfer 
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from the positron to the hadronic system, written 
as a positive quantity; = —q^ = —{k — fc')^, 
where fc, k' are the initial and final state positron 
4-momenta. The Bjorken scaling variable x is 
defined by a: = , where p is the initial proton 
beam 4-momentum. x is the fractional momen- 
tum carried by the struck parton, in the parton 
model of DIS processes. W is the invariant mass 
of the final state hadron system, or the 7*7? invari- 
ant mass given by = {p + ■ These three 
quantities are related by = ^^("'^)- Fi- 
nally, the variable U — ^ — ^ where -y/s is the 
total ep centre-of-mass energy. 

This review is divided into two main sections, 
and a short summary. In section I discuss charm 
photoproduction measurements and QCD calcu- 
lations. Charm production in DIS, and measure- 
ments of the charm contribution to the proton 
structure function, are presented in section ^. 

2. OPEN CHARM IN PHOTOPRODUC- 
TION 

The techniques used by the two collabora- 
tions!^,^ to reconstruct charm candidates in pho- 
toproduction events are similar. The aim is to 
identify the decay chain 

i^*(2010)+ D\t K- 7r+ (1) 

and its charge conjugate, where 7r+ is a "slow" 
pion carrying only 40 MeV of momentum in the 
D*^ rest frame. The analysis uses pairs of op- 
positely charged tracks to search for candi- 
dates, and ZEUS also has a signal from 4-track 
combinations. Each track in turn is assigned the 
charged kaon mass, with the remainder assumed 
to be pions. Combinations whose invariant mass 
rriK-n falls within a window around the mass 
are then further combined with a track of oppo- 
site charge to the kaon candidate, assumed to be 
the slow pion. The number of D*^ decays in 
the sample is estimated from the distribution of 
Am = mK-KTrs ~ f^K-Tr- The mass difference dis- 
tribution from HI is shown in Figure ^ and a 
clear signal can be seen with a relatively small 
background. 

Kinematic cuts are applied to the D* candi- 
dates, to ensure clean reconstruction and reduce 
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Figure 2. Distribution of Am for D* candidates 
in photoproduction, as measured by HI. 



background. Both experiments require their D*s 
to be in the central region of pseudorapidity, 
-1.5 < 7/ < 1.0 (40° < e < 155°). HI apply 
a pt cut of 2.5 GeV, while ZEUS apply cuts of 
3(4) GeV for KtttTs {K-k-k-k-Ks) signals. 

Photoproduction event samples with < 

4 GeV^ and W in the range 100 to 200 GeV are 
identified by activity in the main detector, with 
no scattered positron detected. The average value 
of is around 0.2 GeV^, and the median value 

5 • 10^^ GeV'^. HI additionally use a tagged sam- 
ple where the positron is detected in the small 
angle luminosity detector, giving an average 

of 10-^GeV^ 

From the number of reconstructed _D*s, the 
experiments measure the visible cross-section for 
the reaction e^p — > D* X within the kinematic 
cuts applied, by correcting for acceptance and 
for the combined branching ratio to the decay 
channels detected. In order to extract the total 
charm production cross-section, the visible frac- 
tion falling within the cuts must be estimated 
from Monte Carlo. A further correction must 

be applied for the probability p^ for the 

charm quark to fragment to a D* meson. The 
Monte Carlo extrapolation to the full D* kine- 
matic region is subject to large uncertainties, aris- 
ing mainly from the choice of structure functions 
and of the value of rric used in the calculations, 
and this uncertainty is reflected in the cross- 
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• III tagged data: PREL 

• HI tagged data 

• H 1 untagged data 




NLO (MRSG + GRV) 
upper p/m^ = 0.5 
lower \ijm^ = 2 



Figure 3. The total cross-section for charm pho- 
toproduction, plotted as a function of W. Er- 
ror bars include uncertainties due to the model- 
dependent extrapolation outside the experimen- 
tal acceptance. Data points from HI are plotted 
with inner error bars showing the contribution 
from statistical and experimental systematic er- 
rors only. 



section error bars as plotted in Figure I The 
data show an order of magnitude rise in the inte- 
grated cross-section relative to lower energy pho- 
toproduction measurements, in general agreee- 
ment with the NLO QCD calculations of [Q. 

The calculations of differential cross-sections, 
covering the {pt,r],W) acceptance of the mea- 
surements, are much less sensitive than the to- 
tal cross-section to the uncertainties mentioned 
above. When these distributions are plotted, 
the data are found to lie systematically above 
the predictions of Q. The predictions are based 
on an approach to performing NLO calculations 
where only light quarks are assumed to be active 
flavours in the structure functions of the proton 
and photon. Alternative calculations have been 
performed 1^,^ where charm flavour excitation is 
allowed to contribute in leading order, giving gen- 
erally higher predicted cross-sections. These are 
due principally to the process eg — *■ eg, where 
the hard scattering is between a charm quark 
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ZEUS 94 D -> (K n 71 Jt) 
• ZEUS 94 D* -> (K It) 11^ 

□ HI D* ^ (K 7t) Jl, 




massive 

(tg = m^, m_, - 1.5 GeV 
massive 

[Ig = 0.5 m^, m^, = 1.2 GeV 
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Figure 4. The D* pt spectrum in photoproduc- 
tion compared with QCD calculations. 



from the photon and a gluon from the proton. 
The two types of calculation are described as 
massive and massless charm approaches, respec- 
tively. Figure ^ shows the measured pt spectrum 
from both experiments, together with the massive 
charm predictions of |^ using two different values 
for the charm quark mass mc, and the massless 
charm predictions of Q . The massless charm ap- 
proach is seen to give a better description of the 
data. The calculations of ||^, not plotted in Fig- 
ure |[ are based on a different treatment of the 
heavy quark fragmentation and also agree well. 

3. OPEN CHARM IN DEEP INELASTIC 
SCATTERING 

Measurements of charm meson production in 
DIS events have been used by H10 and Zeus[|j 
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to evaluate the charm contribution to the proton 
structure function F2{x,Q^). This contribution 
F2'^{x, Q^) is defined by 



dxdQ^ 



(l + (f-yf)Ff(x,Q2) (2) 



where a'^'^ is the total charm production cross- 
section. Here the contribution due to longitudi- 
nally polarised virtual photons, which is expected 
to be small, has been neglected. The measure- 
ments of ^2'^ are also discussed in the contribu- 
tion on structure functions to this conference 

The experiments perform their analyses start- 
ing from samples selected by the standard cuts 
for DIS event samples, principally by requiring 
the detection of an energetic scattered positron. 
ZEUS use events in the range 5 < < 100 GeV^ 
and HI require 10 < < 100 GeV^ The D* 
decay channel (|l|) is again used as the signal for 
charm production by both experiments, while HI 
have additionally analysed measurements of in- 
clusive production. 

In order to throw light on the underlying pro- 
duction mechanism, the fractional momentum xb 
of the charm mesons in the 7*p frame is studied. 
This quantity is defined by HI as 



xd 



= 2 



\Pd[ 
W 



(3) 



where and p* are the momenta in this frame 
of the D'^ and the proton, respectively. The dis- 
tribution is shown in Figure ZEUS present a 
similar distribution using the momentum of the 
D* instead of the 13°. If the charm particles are 
produced via photon-gluon fusion (PGF), a cc 
pair is recoiling against the proton remnant in 
the 7*p frame. In this case the xd distribution 
is expected to peak below xd — 0.5, as observed. 
In the flavour excitation process where the vir- 
tual photon interacts with a charm quark from 
the sea, the distribution is expected to be shifted 
towards larger values oi xd- As shown in Fig- 
ure 1^, the data are consistent with expectations 
from a pure PGF Monte Carlo. The distributions 
of the observed events in , rj^ , W and are 
also found to be in good agreement with predic- 
tions. 
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Figure 5. Distribution of D'^ fractional momen- 
tum Xd vo- DIS measured by HI (points) com- 
pared with expectations from a PGF Monte Carlo 
(AROMA[p^, shaded histogram) and from charm 
quark sea contributions (dashed histogram and 



The measurements of ¥2'^ are shown in Fig- 
ure ^, together with results from the EMC fill] 
fixed target experiment at higher x. For the low 
values of x accessible at HERA, x ~ 10~^, it is 
found that i^j'' is around 25% of the total F2. 
Also shown in Figure ^ is the range of predictions 
from an NLO calculation |l2| using the GRV||l|] 
gluon distribution, for charm quark masses be- 
tween 1.35 < rric < 1.7 GeV. The lowest value of 
the mass corresponds to the upper boundary of 
the shaded region in the figure. 

The NLO QCD calculation of [|| includes only 
the PGF diagram as the source of charm produc- 
tion in DIS. This approach is therefore similar to 
the massive charm calculations for photoproduc- 
tion discussed in section ^. There it was seen that 
such an approach gives cross-section predictions 
which tend to be lower than the photoproduction 
data, while for DIS the agreement is good within 
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experimental and theoretical activity. Models 
of photoproduction are in reasonable agreement 
with data at current levels of precision, with some 
suggestion of the need for a contribution due to 
"active charm" in the photon. Charm production 
is found to contribute a substantial fraction of the 
proton structure function at small expected 
given the observed steep rise of the gluon distri- 
bution. Precise measurements of Q^) will 
provide a useful independent constraint on fits to 
parton distributions, and a check of different ap- 
proaches to calculations of heavy quark prodution 
in perturbative QCD. 



Figure 6. Measurements of F^" from 1994 HERA 
data, plotted as a function of a; in 4 bins of Q^. 
Also shown are fixed target measurements from 
EMC. The shaded bands are NLO QCD predic- 
tions for a range of charm quark masses. 



the limited statistical precision of the 1994 data 
samples. 

Again for DIS, there exist alternative treat- 
ments within the framework of perturbative 
QCD, which treat the charm quark as an ac- 
tive parton with its own distribution function 
c(x,Q^). The charm distribution is driven by 
evolution from the gluon density, being zero be- 
low some threshold value of order m^. These 
treatments are expected to give improved re- 
sults for large phenomena. A number of au- 
thors have recently discussed ways to im- 
plement such schemes so that the behaviour at 
~ is essentially given by PGF, to give 
a description of charm production valid for all 
scales. Such analyses predict larger values of F2'^ 
at high than expected for pure PGF, a pre- 
diction which should be testable with the higher 
statistics 1995-6 HERA data now being analysed. 

4. SUMMARY 

Charm production at HERA in real and vir- 
tual photon-proton collisions provides an excel- 
lent laboratory for detailed tests of perturbative 
QCD models. This is an area of considerable 
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